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ABSTRACT 

It is generally assumed that the central galaxy in a dark matter halo, that is, the 
galaxy with the lowest specific potential energy, is also the brightest halo galaxy 
(BHG), and that it resides at rest at the centre of the dark matter potential well. 
This central galaxy paradigm (CGP) is an essential assumption made in various fields 
of astronomical research. In this paper we test the validity of the CGP using a large 
galaxy group catalogue constructed from the Sloan Digital Sky Survey. For each group 
we compute two statistics, 1Z and S, which quantify the offsets of the line-of-sight 
velocities and projected positions of brightest group galaxies relative to the other 
group members. By comparing the cumulative distributions of \1Z\ and |<S| to those 
obtained from detailed mock group catalogues, we rule out the null-hypothesis that the 
CGP is correct. Rather, the data indicate that in a non-zero fraction /bnc(-^0 of all 
haloes of mass M the BHG is not the central galaxy, but instead, a satellite galaxy. In 
particular, we find that /bnc increases from ~ 0.25 in low mass haloes (10 12 /i _1 M Q < 
M < 2 X 10 13 /i _1 M Q ) to - 0.4 in massive haloes (M > 5 x 10 13 /i _1 M ). We show 
that these values of /bnc are uncomfortably high compared to predictions from halo 
occupation statistics and from semi-analytical models of galaxy formation. We end by 
discussing various implications of a non-zero /bncWj with an emphasis on the halo 
masses inferred from satellite kinematics. 

Key words: methods: statistical - galaxies: halos - galaxies: kinematics and dynam- 
ics - dark matter - galaxies: clusters: general 



1 INTRODUCTION 

According to the current paradigm of galaxy formation, all 
galaxies form as a result of gas cooling at the centre of the 
potential well of dark matter haloes. Structures form hierar- 
chically, such that smaller haloes merge to form larger and 
more massive haloes. When a halo and its 'central' galaxy 
is accreted by a larger halo, it becomes a subhalo and its 
galaxy becomes a 'satellite' galaxy. In this paradigm, it is 
assumed that ram-pressure and tidal forces strip satellite 
galaxies of their gas reservoir, causing their star formation 
to be quenched shortly after having been accreted. The cen- 
tral galaxy (i.e., the galaxy with the lowest specific poten- 
tial energy), however, continues to accrete new gas, and is 
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also expected to cannibalize some of its satellites. Conse- 
quently, it is generally assumed that the central galaxy is 
the most luminous, most massive galaxy in a dark matter 
host halo, and that it resides at rest at the centre of the 
halo's potential well. Following van den Bosch et al. (2005; 
hereafter vdB05), we will refer to this as the 'Central Galaxy 
Paradigm' (CGP). 

There are numerous areas of astronomy in which the va- 
lidity of the CGP is an essential assumption, although this is 
rarely enunciated. Examples are various techniques to mea- 
sure halo masses, such as satellite kinematics (e.g., McKay et 
al. 2002; van den Bosch et al. 2004; More et al. 2009), weak 
lensing (e.g., Mandelbaum et al. 2006; Johnston et al. 2007; 
Cacciato et al. 2009; Sheldon et al. 2009b), and strong lens- 
ing (e.g., Kochanek 1995; Cohn et al. 2001; Koopmans & 
Treu 2003; Rusin et al. 2003). In addition, the CGP also fea- 
tures in halo occupation modelling, where assumptions have 
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to be made regarding the distribution of galaxies within dark 
matter haloes in order to compute the galaxy-galaxy corre- 
lation function on small scales (e.g., Scoccimarro et al. 2001; 
Sheth et al. 2001; Yang et al. 2003; Zehavi et al. 2005; Zheng 
et al. 2005; Cooray 2005; van den Bosch et al. 2007; Tinker 
et al. 2008), and in algorithms developed to identify galaxy 
groups and clusters in photometric or spectroscopic redshift 
surveys (e.g., Yang et al. 2005a, 2007; Berlind et al. 2006; 
Koester et al. 2007). 

Whether central galaxies actually comprise a special 
population has been debated for many years. However, re- 
cent analyses of galaxies in groups and clusters (e.g., Wein- 
mann et al. 2006; Skibba et al. 2007; von der Linden et 
al. 2007; van den Bosch et al. 2008; Pasquali et al. 2009, 
2010; Skibba 2009; Hansen et al. 2009) and halo model anal- 
yses of galaxy clustering (e.g., Skibba et al. 2006; Cooray 
2006; van den Bosch et al. 2007; Skibba & Sheth 2009) have 
explicitly shown that central galaxies are indeed a distinct 
population and exhibit different properties (e.g., colour, star 
formation activity, AGN activity, morphology, stellar popu- 
lation properties) than satellite galaxies of the same stellar 
mass, and with different dependencies on the mass of their 
host halo. In addition, these studies have shown that cen- 
tral galaxy properties are strongly correlated with halo mass, 
while those of satellite galaxies only reveal a very weak de- 
pendence on the mass of the halo in which they orbit. How- 
ever, it is important to realize that in virtually all of these 
studies, central galaxies are assumed to the the brightest (or 
most massive) halo galaxies. If this aspect of the CGP is not 
correct for a non-negligible fraction of all haloes, the differ- 
ences between centrals and satellites found by these studies 
will have to be considered lower limits. 

The validity of the CGP has been investigated by a 
number of authors. In particular, several recent, observa- 
tional studies have shown that although most brightest halo 
galaxies (hereafter BHGs) are nearly at rest near the cen- 
troid of their group or cluster, or at the peak of the cluster 
X-ray emission, some of them are not (e.g., Beers & Geller 
1983; Malumuth et al. 1992; Bird 1994; Postman & Lauer 
1995; Zabludoff & Mulchaey 1998; Oegerle & Hill 2001; 
Yoshikawa et al. 2003; Lin & Mohr 2004; von der Linden et 
al. 2007; Bildfell et al. 2008; Hwang & Lee 2008; Sanderson 
et al. 2009; Coziol et al. 2009). These studies have focused on 
either cD galaxies in clusters, or on brightest cluster galax- 
ies (BCGs) in general. It has been argued that most cD 
galaxies form a subpopulation of BCGs (e.g., Bernstein & 
Bhavsar 2001; Coziol et al. 2009) and may have grown by 
"cannibalizing" smaller neighboring galaxies. 

In an early study, Beers & Geller (1983) analyzed the 
spatial distribution of bright galaxies in 56 rich clusters and 
argue that cD galaxies tend to lie at local density peaks but 
not necessarily at the bottom of the potential well of the 
whole cluster. Oegerle & Hill (2001) found that, out of their 
sample of 25 Abell clusters, the cD galaxies of four of them 
have significant peculiar velocities relative to the cluster ve- 
locity. More recent studies have analyzed the positions and 
velocities of BCGs. For example, in a study of 833 SDSS 
clusters, von der Linden et al. (2007) found that 21 BCGs 
in their sample lie further than 1 Mpc away from the mean 
galaxy in the cluster. Hwang & Lee (2008) found that the 
BCGs of two clusters out of a sample of 24 have signifi- 
cantly offset velocities and positions; these two clusters also 



appear to be in dynamical equilibrium. Recently, Coziol et 
al. (2009), using a sample of 452 Abell clusters selected for 
the likely presence of a dominant galaxy, estimated that the 
BCGs have a median peculiar velocity of 32% of their host 
clusters' radial velocity dispersion. 

We emphasize that these results based on large clusters 
do not necessarily hold for less massive haloes. After all, in 
galaxy clusters, the ratio L2/L1 of the luminosities of the 
two brightest galaxies tends to be much smaller in a Milky 
Way (MW)-sized halo, on average (e.g., van den Bosch et al. 
2007). For example, for the halo hosting the MW it is the 
ratio of the luminosities of the Large Magellanic Cloud and 
the MW, which is ~ 0.1 (van den Bergh 1999), while this ra- 
tio is much closer to unity in Virgo, Coma and other nearby 
clusters (Postman & Lauer 1995). Consequently, it is only 
natural to expect that the CGP is more likely to be valid for 
low mass haloes than for massive cluster-sized haloes. Nev- 
ertheless, vdB05 used a large sample of 3473 galaxy groups 
from the group catalogue of Yang et al. (2005a), and con- 
clusively falsified the assumption that the brightest group 
galaxy is always at rest at the centre of the potential well. 
In this paper we extend the analysis of vdB05 using a larger, 
more accurate group catalogue, and focusing on different as- 
pects of the CGP as a function of halo mass. In particular, we 
separately test two aspects of the CGP, namely, (i) central 
galaxies reside at rest at the centre of their host halo's po- 
tential well, and (ii) central galaxies are the brightest, most 
massive galaxies in their host haloes. We do so by comparing 
three hypotheses: 

• Ho'- Our null hypothesis is that the CGP is correct: cen- 
tral galaxies are always the brightest objects in their haloes 
and are at rest at the centre of the potential well. 

• Hi: Central galaxies are the brightest objects in their 
haloes, but they have a velocity and spatial offset with re- 
spect to the centre of the potential well, such that the sys- 
tems still obey the Jeans equations (i.e., they are still in 
dynamical equilibrium). We will specify the amount of off- 
set via a velocity bias parameter, fo vc i, to be defined in Sec- 
tion g3] 

• H2: Central galaxies reside at rest at the centre of the 
potential well, but they are not the brightest objects in a 
fraction /bnc (for 'Brightest-Not-Central') of all dark mat- 
ter haloes. 

In order to avoid confusion, throughout this paper we use the 
term 'central galaxy' to refer to the galaxy with the lowest 
specific potential energy. The central galaxy is the BHG in 
Ho and Hi, but not in H2, and its location coincides with 
the centre of the halo's potential well in Ho and H2 , but not 
in Hi. 

We base our study on the galaxy group catalogue of 
Yang et al. (2007), extracted from the Sloan Digital Sky 
Survey (SDSS; York et al. 2000) Data Release 4 (DR4; 
Adelman-McCarthy et al. 2006). We analyze both the posi- 
tions and velocities of BHGs relative to those of the other 
member galaxies, and compare the results to mock group 
catalogues that correspond to one of our three hypotheses. 
The large SDSS group catalogue, with accurate halo mass 
estimates, allows us to falsify Ho and to quantify the degree 
to which Hi and H2 are valid (i.e., to constrain the values 
of 6 vo i and /bnc). 

This paper is organized as follows. In Section [2] we 
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present two statistics that quantify the offsets of BHGs, 
and which can be used to assess deviations from the central 
galaxy paradigm. We describe the SDSS group catalogue in 
Section [3] and the construction of mock group catalogues in 
Section[4] In Section^ we compare these mock catalogues to 
the data, in order to test our three hypotheses. We find that 
both Ho and Hi can be ruled out, but that Hi yields results 
in agreement with the data as long as 0.25 < /bnc ,$ 0.4, 
with a weak dependence on halo mass. In Section [6] we dis- 
cuss our results and compare them to predictions from halo 
occupation statistics and from two semi-analytic models of 
galaxy formation. We also discuss the implications for stud- 
ies of satellite kinematics. Finally, we end the paper with 
our conclusions and a discussion of additional implications 
of our results. 

Throughout this paper, we adopt a fiat ACDM cosmol- 
ogy with n m = 0.238, Q A = 1 - fim, n = 0.951, a$ = 0.744, 
and we express units that depend on the Hubble constant 
in terms of h = //o/100 kms -1 Mpc -1 . In addition, we use 
'log' as shorthand for the 10-based logarithm. 



2 PHASE-SPACE STATISTICS OF CENTRAL 
GALAXIES 

In order to test the CGP (i.e., hypothesis Ho), we use the 
line-of-sight velocities of galaxies obtained from their red- 
shifts. In what follows, vbhg refers to the line-of-sight ve- 
locity of the BHG, and t) sa t,i refers to the line-of-sight ve- 
locity of the i th satellite galaxy. We define the difference 
AV = Wsat — vbhg between the mean velocity of the satellite 
galaxies and that of the BHG. If the CGP is correct and 
«sat,i follows a Gaussian distribution with velocity disper- 
sion (Tsat, then the probability that a halo with jV sat satellite 
galaxies has a value of AV is given by 



P(AV)dAV = 



27T<7 



exp 



2a 2 



dAV 



(1) 



where a = cr sa t / V JV aat . Therefore, in principle, one could 
define the parameter 



R 



vbhg ) 



(2) 



and test the CGP by checking whether R follows a normal 
distribution with zero mean and unit variance. However, the 
velocity dispersion <7 sa t is generally unknown, and instead we 
must use its unbiased estimator 



- «sat 

\ J Vsat 1 . 

\ 1=1 



(3) 



Following vdB05, we use the following modified parameter 
as an indicator of the offset between BHGs and satellite 



1Z = 



vbhg ) 



0"s a t 



(4) 



This parameter is similar to the relative peculiar velocities 
of brightest cluster galaxies, |Dpec|/er c i ustcr , used in related 
studies (e.g., Malumuth et al. 1992; Coziol et al. 2009). If 
the null-hypothesis of the CGP is correct, TZ should follow a 
Student t-distribution with v = iV sa t — 1 degrees of freedom. 



Note that P v (TZ) approaches a normal distribution with zero 
mean and unit variance in the limit A^at — > 00. 

In practice, although the galaxy group finder (Yang 
et al. 2005a, 2007) has been thoroughly tested with mock 
SDSS catalogues to reliably identify galaxies residing in the 
same dark matter halo, it is not perfect. In particular, be- 
cause of redshift errors and redshift-space distortions, the 
group finder inevitably selects some interlopers (galaxies 
that are not associated with the same halo). In addition, 
the SDSS suffers from various incompleteness effects. If the 
actual BHG is missed or misidentified, TZ will be measured 
with respect to a satellite galaxy, and \1Z\ will tend to be 
overestimated. The presence of interlopers and incomplete- 
ness effects tend to create excessive wings in the 1Z distri- 
bution, and therefore a direct comparison with the Student 
t-distribution cannot be made. To circumvent these prob- 
lems, we compare the ^-distributions obtained from galaxy 
groups identified in the SDSS to those obtained from groups 
identified in mock galaxy redshift surveys (Section|4|, which 
suffer from interlopers and incompleteness to the same ex- 
tent as the real data. 

We also investigate the spatial offsets of BHGs in this 
paper, and to do so we introduce the following parameter, 
analogous to the parameter 7Z, that quantifies the spatial 
separation between BHGs and satellite galaxies, using their 
projected angular separations perpendicular to the line-of- 
sight: 



S = 



r p ,BHG) 



(5) 



where f p . sa t is the mean projected position of satellite galax- 
ies in a group, in terms of the galaxies' mean right ascensions 
and declinations, and 



'>"p,sat 



^ JVsat 

T / (^p,sat.z Vp, s 



(6) 



are expressed in h kpc in 



Both F„ sat — ?>,bhg and a T a1 
the computation of S. 

For the mean and standard deviation u sa t and <7 sa t in 
the TZ parameter (Eqn. 2]), and r p , sa t and ov p , sat in the S 
parameter (Eqn. [5), we have tested that other estimators 
for these quantities, such as the biweight estimator and the 
gapper (Beers et al. 1990), have an insignificant effect on 
our results. 

Note that, with these definitions, the velocity offset TZ 
and spatial offset S are only valid for groups with three 
or more members, and are not defined for galaxy pairs or 
isolated galaxies. 



3 APPLICATION TO THE SDSS 

We first describe the SDSS galaxy group catalogue in Sec- 
tion 13.11 and then the subset of galaxy groups used in our 
analysis in Section f3. 21 



3.1 Galaxy Group Catalogue 

The analysis presented in this paper is based on the SDSS 
galaxy group catalogue of Yang et al. (2007; hereafter Y07), 
which is constructed by applying the halo-based group 



© 2009 RAS, MNRAS 00Q.rTlfl7l 



4 Skibba et al. 



finder of Yang et al. (2005a) to the New York University 
Value-Added Galaxy Catalog (NYU-VAGC; see Blanton et 
al. 2005), which is based on the SDSS DR4 (Adelman- 
McCarthy et al. 2006). From this catalogue Y07 selected 
all galaxies in the Main Galaxy Sample (Strauss et al. 2002) 
with an extinction corrected apparent magnitude brighter 
than m r = 18, with redshifts in the range 0.01 < z < 0.20 
and with a redshift completeness C z > 0.7. 

This sample of galaxies is used to construct three group 
samples: sample I, which only uses the 362356 galaxies with 
measured redshifts from the SDSS, sample II which also in- 
cludes 7091 galaxies with SDSS photometry but with red- 
shifts taken from alternative survey^], and sample III which 
includes an additional 38672 galaxies that lack a redshift 
due to fiber collisions, but which we assign the redshift of 
its nearest neighbor. The present analysis is based on the 
galaxies in sample II with m r < 17.77, which consists of 
344010 galaxies. 

All the magnitudes and colours of the galaxies are Pet- 
rosian, and they have been corrected for Galactic extinc- 
tion (Schlegel, Finkbeiner & Davis 1998), and have been 
fc-corrected and evolution-corrected to z = 0.1 (Blanton & 
Roweis 2007). Stellar masses for all galaxies are computed 
using the relations between stellar mass-to-light ratio and 
colour of Bell et al. (2003). 

The geometry of the SDSS used for the group catalogue 
is defined as the region on the sky that satisfies the redshift 
completeness criterion. To account for the effects of the sur- 
vey edges, Y07 used the SDSS DR4 survey mask with mock 
galaxy redshift surveys to estimate the fraction of "missing" 
members within the halo radius for each group. The group 
luminosities and masses are corrected for this fraction, and 
groups missing 40% or more of their members were excluded, 
which removes only 1.6% of all groups. 

As described in Y07, the majority of the groups in our 
catalogue have two estimates of their dark matter halo mass: 
one based on the ranking of its total characteristic luminos- 
ity, and the other based on the ranking of its total character- 
istic stellar mass, both determined from group galaxies more 
luminous than M r — 5 log h — —19.5. As shown in Y07, both 
halo masses agree very well with each other, with an aver- 
age scatter that decreases from ~ 0.1 dex at the low mass 
end to ~ 0.05 at the massive end. In this paper we adopt 
the group masses based on the stellar mass ranking, but we 
have checked that the luminosity ranking gives results that 
are almost indistinguishable. The stellar mass based group 
masses are available for a total of 215493 groups in our sam- 
ple, which host a total of 277838 galaxies. This implies that 
a total of 66172 galaxies have been assigned to a group for 
which no reliable mass estimate is available (but see Yang, 
Mo & van den Bosch 2009a). 



bers, with 50 km s _1 < <r sa t < 1000 kms -1 , and with reliable 
group masses greater than 10 12 ft -1 Mq. 

Because of the finite thickness of the spectroscopic fibers 
used, the SDSS suffers from incompleteness due to fiber col- 
lisions. No two fibers on the same SDSS plate can be closer 
than 55 arcsec. Although this fiber collision constraint is 
partially alleviated by the fact that neighboring plates have 
overlap regions, ~ 7 percent of all galaxies eligible for spec- 
troscopy do not have a measured redshift. Since fiber colli- 
sions are more frequent in regions of high (projected) den- 
sity, they are more likely to occur in richer groups, thus 
causing a systematic bias that may need to be accounted 
for. Although Sample III tries to correct for this incom- 
pleteness by assigning galaxies that lack a redshift due to 
fiber collisions the redshift of its nearest neighbor, Zehavi et 
al. (2002) have shown that in roughly 40 percent of cases, 
the redshift thus assigned carries a large error. Hence, al- 
though Sample II is more incomplete than Sample III, there 
is less of a risk of interlopers, and the redshifts are accurate, 
which is necessary for the analysis with galaxy line-of-sight 
velocities. However, if the true brightest galaxy in a group is 
missed due to a fiber collision, the velocity and spatial off- 
sets 1Z and S will be estimated relative to a satellite galaxy, 
and will tend to be overestimated, resulting in a stronger 
signal. To avoid this problem, we exclude groups in Sample 
II in which the brightest galaxy is not also a brightest group 
galaxy in Sample III. 

This results in a sample of 6760 groups with iV ga i > 
3, excluding approximately 7% of the galaxy groups. This 
constitutes our fiducial galaxy group sample0 To be clear, 
our fiducial sample simply consists of a subset of Sample II 
groups. We emphasize that Sample III groups are not used in 
our analysis; they are only used to remove those groups from 
Sample II that may have been affected by fiber collisions. 

In Fig. [T] we show the abundance of galaxy groups as a 
function of richness in the catalogue, with M > 10 12 /i _1 M© 
and 50 kms -1 < frrmsat < 1000 kms -1 . The requirement 
that the brightest group galaxy is also the brightest group 
galaxy in Sample III results in slightly fewer groups at all 
richnesses. 

The analysis of vdB05 was done with a catalogue of 2502 
groups in the 2dFGRS with four members or more; the new 
Y07 catalogue has nearly twice as many groups with 7V ga i > 
4 (4571), which is a significant improvement. In addition, the 
typical rms redshift and magnitude errors of galaxies in the 
SDSS are 30 kms -1 and 0.035 mag (r-band), respectively 
(Strauss et al. 2002), compared to 85 kms -1 and 0.15 mag 
(B-band) in the 2dFGRS (Colless et al. 2001). The improved 
statistics of this SDSS DR4 group catalogue allow us to not 
only investigate the halo mass dependence of central galaxy 
velocity bias, but also the dependence of velocity bias on the 
properties of the central galaxies themselves. 



3.2 Galaxy Groups used in this Paper 

In what follows we restrict our analyses to the 7234 galaxy 
groups in the sample II catalogue with three or more mem- 



1 These redshifts are taken from the 2dFGRS (Colless et al. 
2001), IRAS PSCz (Saunders et al. 2000), or RC3 (de Vau- 
couleurs et al. 1991). See Blanton et al. (2005) for details. 



2 Using the additional criterion that the group masses in Samples 
II and III are within 0.3 dex of each other excluded another 7%, 
but yielded results that were indistinguishable from those based 
on our fiducial sample. 
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4.1 Assigning Luminosities 



1000 



100 




Figure 1. Galaxy group multiplicity function of SDSS Sample II. 
Circle points show the number of groups as a function of group 
richness, for groups with M > 10 12 h~ 1 Mq and 50 kms -1 < 
CTsat < 1000 kms -1 . Square points show JV gr p(iVg a i) for groups 
in which we additionally require that the central galaxy is also 
the central galaxy of a group in Sample III. Horizontal error bars 
indicate the width of each bin. 



4 MOCK CATALOGUES 

The main goal of this paper is to use the distributions of 
the parameters 1Z and <S defined above to test our three 
hypotheses related to the CGP defined in Section [T] As dis- 
cussed above, since the SDSS galaxy group catalogue suffers 
from interlopers and incompleteness effects, we require mock 
group catalogues constructed from mock galaxy redshift sur- 
veys (hereafter MGRSs) using the same halo-based galaxy 
group finder as used for the SDSS. 

We construct MGRSs by populating dark matter haloes 
with galaxies of different luminosities. The distribution of 
dark matter haloes is obtained from a set of large iV-body 
simulations (dark matter only) for the WMAP3 ACDM cos- 
mology from Maccio et al. (2007). The simulations have 512 3 
particles each, have periodic boundary conditions, and box 
sizes of Lbox = 100/i -1 Mpc (hereafter Lioo) and I/box = 
300/T 1 Mpc (hereafter L 3 oo)- We follow Yang et al. (2004) 
and replicate the L300 box on a 4 x 4 x 4 grid. The central 
2x2x2 boxes, are replaced by a stack of 6 x 6 x 6 L100 boxes 
(see Fig. 11 in Yang et al. 2004). This stacking geometry cir- 
cumvents incompleteness problems in the mock survey due 
to insufficient mass resolution of the L300 simulations, and 
allows us to reach the desired depth of z ma x = 0.20 in all 
directions. 

Dark matter haloes are identified using the standard 
FOF algorithm with a linking length of 0.2 times the mean 
inter-particle separation. Unbound haloes and haloes with 
less than 10 particles are removed from the sample. The 
resulting halo mass functions are in excellent agreement with 
the analytical halo mass function of Sheth, Mo & Tormen 
(2001). 



We populate each halo with galaxies of different luminosities 
using the conditional luminosity function (CLF) model de- 
scribed in Cacciato et al. (2009; hereafter C09). The CLF, 
&(L\M), specifies the average number of galaxies of lumi- 
nosity L in a halo of mass M, and is constrained to accu- 
rately match the SDSS r-band luminosity function (Blanton 
et al. 2003), the clustering strength of SDSS galaxies as a 
function of luminosity (Wang et al. 2007), and the galaxy- 
galaxy lensing data of Mandelbaum et al. (2006). The CLF 
assumes that the luminosity-dependent abundance, distribu- 
tion, and clustering of galaxies can be described as a func- 
tion of halo mass. The CLF of C09 is split in two parts, 
$cen(L|M) and $ sa t(LjA/), which describe the halo occupa- 
tion statistics of central and satellite galaxies, respectively. 

For each halo we draw the luminosity of its central 
galaxy from <l> CGn (L|M), which is parameterized as a log- 
normal distribution: 



<!>con(L|M)dL 



27rln(10)cr c . 



■ exp 



\/2cr ccr 



AL 
T 



(7) 

Here cr ccn = 0.14 quantifies the scatter between central 
galaxy luminosity and host halo mass. More et al. (2009) 
obtained a similar value from their analysis of satellite kine- 
matics: (j ccn = 0.16 ± 0.04. The central galaxy luminosity 
as a function of halo mass, L cen (M) is parameterized as a 
double power-law, with a slope of 3.3 in low-mass haloes and 
0.26 in massive haloes (see C09 for details). 

For the satellite galaxies we assume that their halo oc- 
cupation numbers follow a Poisson distribution with mean 



(N sat \M) 



$ sat (L|M) dL, 



(8) 



where we adopt a luminosity threshold, L m i n , corresponding 
to M r — 51ogh, = — 14. The satellite luminosities are drawn 
from the satellite CLF $ sat (L|M), which is parameterized 
as a modified Schechter function: 



$ sat (L|M) = ^ — ] exp 



L 



(9) 



where </> sat and a sa t are functions of halo mass, the param- 
eter s = 2, and Lj at (M) = 0.562 L ccn (M) (see C09 for de- 
tails). We emphasize that this particular form for the CLF 
(eqn. [7] and [9]) is not an assumption, but rather is the form 
that agrees with the CLF obtained directly from the SDSS 
group catalogue (see Yang et al. 2008). The only poorly 
constrained parameter is s, and we test the effect of its un- 
certainty in Section T6. II 

If the luminosity of the satellite is brighter than that 
of its central, a new luminosity is drawn until it is fainter 
than that of the central. Hence, in our mock universe central 
galaxies are always BHGs, by construction. 



4.2 Assigning Phase-Space Coordinates 

Having assigned all mock galaxies their luminosities, the 
next step is to assign them a position and velocity within 
their halo. 

We assume that each dark matter halo of mass M has a 
NFW (Navarro, Frenk & White 1997) density distribution, 
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Pdm(r\M), with virial radius r v i T (M), characteristic scale ra- 
dius r s (M), and concentration parameter c(M) = r v i t /r s . 
We model the halo concentrations using the c(M) relation 
of Maccio et al. (2007). Assuming haloes to be spherical 
and isotropic, the local, one-dimensional velocity dispersion 
follows from solving the Jeans equation 



' p dm (r'\M)^(r'\M)dr' (10) 



0dm(r\M) = . . . 

Pdm(r|M) 

with \?(r) the gravitational potential (Binney & Tremaine 
1987). Using that d^ /dr = G M(r)/r 2 and defining the 
virial velocity Vvit = y/ GM/r v - 1T we obtain 

aL(r|M) = VX-^r (-) ( 1 + - V Ar/rs) (11) 



(12) 



" f( c ) \ r s ) \ r s 
with f(x) = ln(l + x) — x/(l + x) and 

/(r)dr 



,„ r3(l + r)2 

The halo-averaged velocity dispersion is given by 

rr vir (M) 



<a dm |M> = — 
= Ki, 



Pdm(r\M) G dm {r\M)r' dr 

y s/2 z 1/2 (y) 



(i + v) 



■dy 



(13) 



/ 3 (c) ,0 
(cf. van den Bosch et al. 2004). 

4-2.1 Central Galaxies 

For the central galaxies, we proceed as follows. In the case 
of Ho and Hi mocks, we position the central galaxy at rest 
at the centre of the dark matter halo (for the Hi mocks, we 
then reshuffle the indices of central and brightest satellite, 
as detailed below, but we do so only after the construction 
of the group catalogue). In the case of the Hi mocks, we 
follow the approach of vdB05, to which we refer the reader 
for details. Briefly, we assume that the radial coordinate of 
the central, r, follows a probability distribution^ 

2 



Pcen(r)dr : 



(r + a 



^dr, 



(14) 



where a is a free parameter, which is related to the veloc- 
ity bias fovoi as detailed below. In order to parameterize the 
characteristic radius a in terms of that of the dark matter 
halo, we define the parameter / con = a/r B . 

Central galaxies in a halo of mass M at a halo-centric 
radius r have an isotropic velocity dispersion 

1 



,(r|M) = 



,(r|M) 



' P cen(r'\M)^-(r'\M)dr' 



dr 



V 2 



f(c) 



fcen ~h 



with 



J{y) 



fir) dr 

r 3 (/can+T) 



J(r/r s ) (15) 



(16) 



3 The choice for this particular probability distribution, which 
corresponds to a Hernquist (1990) profile, is not motivated by 
any physical considerations, other than the fact that it is well 
behaved, both at r = and at r — > 00. Our results do not depend 
significantly on the exact shape of this probability distribution. 
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Figure 2. The velocity bias (solid curve) and spatial bias (dashed 
curve) of central galaxies as a function of the parameter /cen, 
which expresses the characteristic scale of the radial distribution 
of central galaxies in terms of the characteristic scale of the NFW 
density distribution. The results shown correspond to a dark mat- 
ter halo with a concentration c = 10. 



This implies a halo-averaged velocity dispersion of 



{a cen \M} = 



J r ™ {M> Pcen{r\M) <Xccn(r|M) r 2 dr 
J ^ iAM) Pccn (r\M)r* dr 

y 3/2 J 1/2 (y) 



Kir 



Ac 



dy, (17) 



f(c) J ^Jo (/ccn + y) 3 / 2 

which allows us to define the velocity bias of central galaxies 



en\M) 



en|M) 



(18) 



(a dm \M) (a sat |M) • 
In addition to the velocity bias, we define the spatial bias as 
, _ (r oen |M) (r cen |M) 



" <r dm jAf> (r sat |M) ' 
where the expectation value for the radius follows from 

f rvir p(r)r 3 dr 



(19) 



(r\M) = 



JT™ p(r)r 2 dr 



(20) 



For comparison, an NFW density distribution has (r\M) = 
0.41r V ir(M) for c = 10, and 0.47r vir (M) for c = 5. 

With this model, a particular value of the parameter 
/ ccn implies a particular amount of velocity and spatial bias. 
The relations 6 V ei(/cen) and 6 ra d(/cen) are shown in Figure^ 
These relations only depend weakly on halo concentration 
(see vdB05). In the limit / CGn — > 0, the probability distri- 
bution P en(r) becomes a delta function, implying that the 
central galaxy is sitting at rest at the centre of the dark mat- 
ter halo (i.e., the null-hypothesis Ho of the CGP). Larger 
values of / con result in larger amounts of velocity and spa- 
tial bias. Note that b VB \ is always larger than 6 ra d, indicating 
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that the signature of an off-centered central galaxy (i.e., hy- 
pothesis Hi) is more pronounced, and thus easier to detect, 
in velocity space than in configuration space. 

4-2.2 Satellite Galaxies 

Throughout this paper, we assume that the iV sat satellite 
galaxies in a halo of mass M follow a number density dis- 
tribution n sa t(f|M) = (N sat /M)p dm (r\M), so that there is 
no spatial bias between satellite galaxies and dark mat- 
ter particles. If we further assume that the satellites are 
in isotropic equilibrium, it also follows that there is no ve- 
locity bias between the satellites and the dark matter, nei- 
ther globally [i.e., (cr S at|M) = {adm\M}] nor locally [i.e. 
cr sat (r|M) = a dm (r|M)]. 

For simplicity, we assume that the satellites follow a 
spherically symmetric spatial distribution with a velocity 
distribution that is locally isotropic. Although this is a 
clear oversimplification, as some galaxy groups and clus- 
ters have anisotropic or aspherical distributions (e.g., Bailin 
et al. 2008; Wang et al. 2008), we do not believe that it 
strongly impacts our results. For example, the global ve- 
locity dispersion (i.e., obtained from all satellites) of an 
anisotropic system will be nearly identical to that of an 
isotropic system with the same gravitational potential, since 
both are governed by the virial equation. In other words, 
anisotropy changes the local line-of-sight velocity distribu- 
tion (LOSVD), but leaves the second moment of the global 
LOSVD largely unchanged (see, e.g., van den Bosch et 
al. 2004). Therefore, since the TZ parameter (eqn.[4]| only de- 
pends on the first and second moments of the global LOSVD, 
our results are robust to the assumptions of isotropy and 
sphericity. The cumulative TZ distribution is weakly depen- 
dent on higher moments of the LOSVD, but again, the effect 
of anisotropy is very small. 

Although there is evidence to support our assumption 
that satellite galaxies have number density distributions 
that are well fit by a NFW profile (e.g., Lin et al. 2004), 
the corresponding concentration parameters seem to be sig- 
nificantly smaller than the values expected for their dark 
matter haloes, by about a factor of 2 to 3 (e.g., Hansen et 
al. 2005; Yang et al. 2005b). Since we assume that satellites 
are unbiased with respect to the dark matter, our value for 
6vci will be incorrect by a factor of (<T sa t)/(cdm). Using the 
model described in Eqs (|13[) and (|17[1 , we estimate that this 
may result in an underestimating of the velocity bias by no 
more than 5%. Given that this is a small effect compared to 
the other uncertainties involved in our analyses, we do not 
attempt to correct for it. 

4.3 Mock Group Catalogues 

Once the dark matter haloes are populated with galaxies, 
we construct a MGRS following the procedure described in 
Li et al. (2007). We place a virtual observer at the cen- 
tre of the stack of simulation boxes, assign each galaxy a 
(a, <5)-coordinate, and remove the ones that are outside the 
mocked SDSS survey region. For each model galaxy in the 
survey region, we compute its redshift (which includes the 
cosmological redshift due to the universal expansion, the 
peculiar velocity, and a 35kms _1 Gaussian line-of-sight ve- 
locity dispersion to mimic the redshift errors in the data), 



and its r-band apparent magnitude (based on the r-band 
luminosity of the galaxy). 

We eliminate galaxies that are fainter than the SDSS 
apparent magnitude limit, and incorporate the position- 
dependent incompleteness by randomly eliminating galax- 
ies according to the completeness factors obtained from 
the survey masks provided by the NYU-VAGC (Blanton et 
al. 2005). Finally, we construct group catalogues from the 
MGRSs, using the same halo-based group finder as used for 
the real SDSS DR4. 

Using the method outlined above, we construct the fol- 
lowing set of mock group catalogues. The first is a set of 
ten mocks that only differ in the value of the velocity bias 
6 vc i = 0.0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8, and 1.0. These 
constitute our Hi mocks. Note that the mock with b va \ = 
satisfies the CGP, and therefore corresponds to the null- 
hypothesis Ho, while the other mocks correspond to Hi (i.e., 
BHG is central galaxy, but is not at rest at centre of dark 
matter potential well). The second set of mock group cata- 
logues is constructed starting from the 6 vc i = mock group 
catalogue, in which we switch the luminosities of the cen- 
tral and its brightest satellite for a random fraction /bnc 
of all groups. We construct 6 mock group catalogues for 
/bnc = 0.1,0.2,0.3,0.4,0.5, and 0.6, respectively. These 
constitute our H2 mocks. 

In order to facilitate the comparison with the SDSS 
galaxy group catalogue, in what follows we discard all mock 
groups with an assigned group mass less than 10 12 /i _1 M©, 
and with a satellite velocity dispersion cr 8 at < 50 kms -1 or 
3"sat > 1000 km s -1 . For the mock groups that have not been 
discarded, we compute the parameters TZ and 5. In the next 
section, we compare the resulting distributions P(< \1Z\) 
and P(< \S\) with those obtained from the SDSS group 
catalogue in order to test our three hypotheses Ho, Hi, and 
Hi. 

Finally, we note that the halo centres and central and 
satellite galaxies are defined only in the mock group cata- 
logues, not in the SDSS group catalogue. The 1Z and 5 pa- 
rameters used in the following analysis are defined with re- 
spect to the BHGs in the mock and observed galaxy groups. 
This way, we are able to statistically analyze groups of a 
wide range in mass, including poor groups for which a cen- 
tre or central galaxy might not be well-defined. 



5 RESULTS 

Fig.[3]shows the cumulative distributions of \1Z\ (left panels) 
and |5| (right panels) obtained from the Hi mocks (upper 
panels) and the Hi mocks (lower panels) for groups with 
four or more members (i.e., iV sa t > 3). As expected, P{\1Z\) 
and P(|5|) become broader for larger values of 6 ve i (upper 
panels) and /bnc (lower panels). A few trends are appar- 
ent. Firstly, note that the |5|-distributions are all bunched 
together for 6 ve i < 0.5. Only for 6 vc i > 0.5 are the cumu- 
lative 1 5 1 -distributions notably different. This shows that it 
is difficult to constrain 6 vc i using the angular positions of 
galaxies, on which the 5-statistic is based, vis-a-vis using 
the velocity-statistic 7Z. Physically this is a reflection of the 
fact that dark matter haloes have steep potential wells, so 
that large velocities are required even for relatively modest 
excursions from the centre. Secondly, and most importantly, 
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Figure 3. The cumulative distributions of \1V\ (left panels) and 
|5| (right panels) for the groups in mock group catalogues with 
four members or more. Upper panels show the distributions for 
mock catalogues with different amounts of b ve i: the solid curves 
are for f> ve i = 0, 0.1, 0.2, 0.3, and 0.4 (red, green, magenta, 
cyan, and blue curves, respectively) and the dashed curves are 
for 6 ve i = 0.5, 0.6, 0.7, 0.8, and 1 (red, green, magenta, cyan, and 
blue curves, respectively). Lower panels show the distributions for 
mock catalogues with different fractions /bnc = 0, 0.1, 0.2, 0.3, 
0.4, 0.5, and 0.6 (red, green, magenta, cyan, blue, dashed red, and 
dashed green curves, respectively). 



Halo Mass Bins 



log M a 



log M log M n 



roup 



low mass 12 12.86 13.3 1434 

intermediate mass 13.3 13.53 13.75 1540 

high mass 13.75 14.09 15.2 1555 

Table 1. Halo mass bins (log(M/h _1 M )) in the SDSS group 
catalogue, for groups with four members or more and 50 km s _1 < 
<5sat < 1000 kms -1 . The log mass ranges, mean masses, and 
number of groups are given. 



comparing the upper and lower panels, it is clear that a non- 
zero value of /bnc has a different impact on P(< \1Z\) and 
P(< \S\) than a non-zero value of b vc \. This is good news, as 
it implies that we will be able to discriminate between our 
three hypotheses. 

In what follows we will compare the cumulative \1Z\ and 
\S\ distributions obtained from the SDSS with those ob- 
tained from our mock catalogues. Since we have a relatively 
large number of groups in our SDSS sample, we can per- 
form this test for a (small) number of bins in group mass, 
thus giving some leverage on a possible halo mass depen- 
dence. At the massive end of the group mass distribution, 
the group catalogue is roughly complete, and the mass distri- 
bution closely follows the halo mass function. The constraint 
of four or more group members, however, cuts off the dis- 
tribution at the low mass end, leaving virtually no groups 



with M < 10 12 /i _1 Mq. We split the SDSS group catalogue 
in three mass bins (listed in Table [TJ that contain a simi- 
lar number (~ 1500) of groups, and the corresponding mass 
cuts occur at fog(M/h _1 Mq) = 13.3 and 13.75. 

The results presented in Sections 15.11 and 15.21 are based 
on groups with at least four members (i.e., with iV S at > 3). 
We have repeated these analyses using samples of groups 
with three, five, six, and ten or more members. We have also 
tested lower and higher halo mass thresholds. In each case we 
obtain constraints on 6 ve i and /bnc that are consistent with 
each other at the la level. Hence, our results do not depend 
significantly on the multiplicity and halo mass thresholds 
used. 



5.1 Testing Hypothesis Hi 

In order to test hypothesis Hi, we compare the SDSS group 
catalogue to mock group catalogues with different values of 
fevci while setting /bnc = 0. Figure [4^, shows the cumulative 
7?.|-distribution for different bins of halo mass obtained from 
the SDSS group catalogue (solid lines) and from three mock 
group catalogues corresponding to different values of 6 V ei 
(0, 0.5, and 1). In all cases only groups with four members 
of more are considered. The numbers in square brackets in 
each panel indicate the range of log(M/h~ 1 Mq) considered. 
Note that at fixed 1Z, the corresponding P(< \1Z\) becomes 
smaller when using a bin with larger group masses. This does 
not necessarily imply a trend of b ve i with halo mass, though. 
It may also be due to a mass dependence of the fraction 
of interlopers, or a mass dependence of the completeness 
of group members. Indeed, the mock group samples for a 
fixed value of fc vo i reveal the same trend, indicating that it 
is most likely an artifact introduced by the group finder. 
For each halo mass bin, the mock catalogue with b vc i — 

(i.e., the one that fulfills the null- hypothesis, Ho, of the 
CGP) predicts a \1Z\ -distribution that is much narrower than 
that of the SDSS, while the mock catalogue with 6 vo i 

1 yields a distribution that is too broad. Interestingly, the 
intermediate case, with b vc i = 0.5, results in P(< \1Z\) that 
are very similar to those of the SDSS, for each mass bin. 
This rules out the null hypothesis Ho, and seems to suggest 
that instead central galaxies have a relatively large velocity 
bias with (cr cen |Af) ~ 0.5(cr sa t|M) with little dependence on 
halo mass M. 

However, Fig.|4p, which is similar to Fig. [4^ except that 
it shows the cumulative |S (-distributions, paints a different 
picture. Both the fe vc i = 0.0 and & vc i = 0.5 mocks predict |5|- 
distributions that are clearly too narrow, for each mass bin. 
Rather, the SDSS |5|-distribution seems to require 0.5 < 
6 vc i < 1.0. Hence, it appears that no single value of b vc \ can 
simultaneously match P{< and P(< \S\) obtained from 
the SDSS, ruling against our hypothesis Hi. In other words, 
although the BHGs have non-zero velocities, the velocity 
spread is too small to explain their displacements from the 
halo centres. 

In order to make this more quantitative, we use the 
Kolmogorov-Smirnov (hereafter KS) test to compute the 
probabilities P K s that P(\Tl\) and P(\S\) of the SDSS and a 
particular mock are drawn from the same distribution. The 
results are shown in Fig. [5] which plots log(PKs) as function 
& ve i for each of the four mass bins considered. Using the me- 
dian and the 16 and 84 percentiles of Pks(6 V oi) for the \1Z\ 
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Figure 4. The cumulative distribution of (left figure) and \S\ (right figure) obtained from the SDSS group catalog (solid black curve) 
compared with those obtained from three of our mocks, with fe vc i = (red dotted curve), 0.5 (blue short-dashed curve), and 1 (green 
long-dashed curve), for groups with four members or more. Results arc shown for four log halo mass intervals, as indicated in square 
brackets in each panel. 




Figure 5. The KS-probability that the cumulative 1Z distribution 
(solid lines) and S distribution (dashed lines) obtained from the 
SDSS groups is consistent with that obtained from our mocks, as a 
function of b ve i. Results are shown for four log halo mass intervals, 
indicated in square brackets in each panel. The horizontal dotted 
line in each panel indicates Pks = 0.01: based on estimates of the 
scatter due to cosmic variance, we consider two distributions to 
be statistically equivalent when Pks > 0.01. 



distributions (solid lines) we obtain fo vcl = 0.47to: 7 

for the 

entire mass range, and 6 vo i = 0.441q or> ^vci ~ 0.50 =t 0.05, 
and fo vc i = 0.43lQ Qg for the low-mass, intermediate-mass, 
and high-mass intervals, respectively. For the |<S| distribu- 
tions (dashed lines) this yields b vc i = 0.82^q'q| for the low- 
mass interval, and fe vc i = 0.831q'q6 f° r the intermediate- and 
high-mass intervals. Therefore, the amounts of velocity bias 
consistent with the 1Z and 5 distributions of SDSS groups 
are significantly different, with their best-fit values of 6 ve i 
differing from each other by more than 4cr. We conclude 
that hypothesis Hi is ruled out, since it cannot simultane- 
ously explain the \1Z\ and \S\ distributions obtained from 
the SDSS group catalogue. 

5.2 Testing Hypothesis Hi 

Having ruled out both Ho and Hi, we now turn our attention 
to hypothesis Hi and investigate whether there is a value of 
/bnc for which the mock group catalogue yields \1Z\- and 
S [-distributions that are consistent with the SDSS data. 
To that extent we compare the SDSS group catalogue to 
mock catalogues with 6 ve i = 0, but with different values of 
/bnc- We proceed in the same way as for Hi in the previous 
section: for each mock, which corresponds to a different value 
of /bnc, we compute P(< \1Z\) and P(< |5|), which we 
compare to the corresponding distributions obtained from 
the SDSS, resulting in a value for Pks, the KS-probability 
that the mock and SDSS distributions are consistent with 
each other. 

The results are shown in Fig. [6] where the solid and 
dashed lines once again correspond to the 1Z and 5 distribu- 
tions, respectively. These KS probabilities are the means of 
log(PKs) for 100 realizations with different random seeds. 
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Figure 6. The KS-probability that the cumulative 72. distribution 
(left figure) and S distribution (right figure) obtained from the 
SDSS groups with four or more members is consistent with that 
obtained from MGRSs, as a function of /bnCi the fraction of 
groups in which the most luminous satellite is brighter than the 
central galaxy (see text for details) . The KS probabilities here are 
the means of log(PKs) for 100 realizations with different random 
seeds. 



The solid lines show that the 1Z distribution of SDSS 
groups indicates a relatively large fraction of groups in 
which the central galaxy is not the BHG, with /bnc in- 
creasing from ~ 0.25 in low mass haloes (10 12 /i -1 Mq < 
M < 2 x 10 13 /i _1 M Q ) to ~ 0.4 in massive haloes [M > 5 x 
10 i3. 75^-1 u[ Q y Interestingly, the KS probabilities obtained 
from the 5 distributions shown by the dashed lines yield 
best-fit values of /bnc that are very similar. Hence, we con- 
clude that contrary to hypothesis Hi, hypothesis H2 can 
simultaneously match the P(< \1Z\) and P{< \S\) obtained 
from the SDSS group catalogue. 

Fig. [7] summarizes our results. It shows the best-fit val- 
ues of 6 vc i (upper panel) and /bnc (lower panel) as func- 
tions of halo mass, as inferred from the |7?.|-distributions 
(squares) and |<Sj-distributions (triangles) obtained from the 
SDSS galaxy group catalogue. The figure clearly shows that 
H2 (i.e., a non-zero /bnc) can simultaneously explain both 
the velocities and the positions of BHGs with respect to 
their satellites, while Hi (i.e., a non-zero b V ei) is ruled out 
because it cannot. 

Considering the \1Z\- and |5|-distributions simultane- 
ously, we infer best-fit values for the fraction of haloes 
in which centrals are not BHGs of /bnc = 26±^% for 
12 < log(M//i _1 M Q ) < 13.3, increasing to /bnc = 30± 3 % 
for haloes with 13.3 < log(M//l _1 M©) < 13.75 and /bnc = 
43^3% in the most massive haloes with log(M/h~ 1 Mq) > 
13.75. As in Section 15.21 the best-fit values and uncertain- 
ties are determined from the median and the 16 and 84 per- 
centiles of the -Pks(/bnc) distributions. 

Before proceeding with a discussion regarding the im- 



Figure 7. Halo mass dependence of b ve \ (upper panel) and /bnc 
(lower panel), inferred from the |72|-distributions (squares) and 
|S|-distributions (triangles) obtained from the analyses of the 
SDSS group catalogue (see text for details). For the results in 
the upper panel, the mock catalogues had different values of fe vo i 
and /bnc = 0; for the lower panel, the mocks had different values 
of /bnc and b vc \ — 0. Points show best-fit values of the parame- 
ters estimated from the KS probabilities, vertical error bars show 
1-cr uncertainties from the 16 and 84 percentiles of the Pks dis- 
tributions, and horizontal error bars indicate the widths of the 
mass bins. Hypothesis Hi is clearly ruled out from the fact that 
the 72-data data implies a velocity bias, 6 ve i, that is inconsistent 
with the value inferred from the 5-data. 



plications of these findings, we briefly address the possibility 
that both Hi and Hi are true; BHGs are not central galax- 
ies in a non-zero fraction /bnc of all haloes, and central 
galaxies have a non-zero value for their velocity bias 6 ve i- 
Using mock group catalogues with non-zero values for both 
/bnc and 6 vc i we estimate that the data can accommodate 
a small amount of velocity bias 6 vc i ,$,0.2. We emphasize, 
though, that the data do not require a non-zero b ve \. 



5.3 The Impact of Substructure 

In the tests described above, we have always assumed that 
satellite galaxies follow a smooth, spherically symmetric, 
number density distribution, and we assigned the satellites 
peculiar velocities within its host halo under the assumption 
that the corresponding potential is smooth. This ignores the 
fact that dark matter haloes are believed to have a wealth of 
substructure (see Giocoli et al. 2010, and references therein; 
on substructure in galaxy clusters, see Richard et al. 2010). 
Satellite galaxies are believed to be associated with these 
substructures. Our treatment of the spatial and kinematic 
properties of satellite galaxies is only consistent with this 
concept of substructure if: (i) each subhalo hosts at most 
one satellite, and (ii) the positions and velocities of sub- 
haloes are not correlated with those of other subhaloes. In 
reality, neither of these criteria is likely to be met. Massive 
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subhaloes are likely to host multiple satellites, and the large 
scale filamentary structure is believed to introduce (weakly) 
correlated directions of infall for subhaloes (e.g., Vitvitska et 
al. 2002; Aubert, Pichon & Colombi 2004; White et al. 2010). 
Both of these effects result in correlations among the posi- 
tions and/or velocities of different satellite galaxies within 
the same host halo, which is likely to have an impact on the 
TZ and S statistics. 

In order to have a crude estimate of the impact of 
substructure, we proceed as follows. We start by popu- 
lating 1000 (virtual) dark matter haloes, all with M — 
3 x 10 14 /i _1 Mq, with galaxies using the same CLF model 
as in Section [3] Similar to the Ho mocks, we assure that 
the central galaxy is always the brightest galaxy in the halo, 
and we position it at rest at the centre of the halo. We 
then compute 5 and TZ for each of these haloes. The black, 
solid lines in Fig. [5] indicate the corresponding P(< \1Z\) 
and P(< \S\). Next we repeat this exercise, but now, in 
each halo, a fraction / su b of all satellites is 'clumped' to- 
gether in a single substructure. We model this substructure 
as a halo of mass m — / su bM (i.e., we assign the galax- 
ies in this subhalo phase-space coordinates in exactly the 
same way as we would if the halo was a 'host' halo of the 
same mass). The phase-space coordinates of the centre of 
the subhalo are drawn in the same way as the phase-space 
coordinates of the satellites that are not in a substructure. 
The long-dashed green curves, short-dashed blue curves, and 
dotted red curves in Fig.[S]show the P(< \TZ\) and P(< \S\) 
thus obtained for / su b = 1/3, 1/5, and 1/10, respectively. 
A comparison with the no-substructure case (black, solid 
line) shows that substructure significantly affects the TZ and 
5 distributions, and hence our conclusions, if all dark mat- 
ter haloes have a most massive substructure whose mass is 
m S ub ^ 0.1M. In particular, if many haloes contain substan- 
tial substructure, then an alternative explanation for the fact 
that our fiducial model (which follows the CGP) does not 
match the data may be that the phase-space coordinates 
of satellite galaxies within the same dark matter host halo 
are correlated. In that case the fractions /bnc obtained in 
Section [5.21 will be significantly overestimated. 

However, using the subhalo mass functions of Giocoli 
et al. (2010), we estimate that only ~ 8% of host haloes 
with M = 3 X 10 14 /i _1 Mq have a most massive subhalo 
with / su b = m/M > 0.1, while only ~ 0.7 percent have a 
most massive subhalo with m/M > 0.3. Based on these es- 
timates, and on the tests described above, we argue that ig- 
noring substructure when populating haloes with (satellite) 
galaxies does not have a significant impact on the TZ and 
S statistics, though more sophisticated tests are required 
to confirm this. Therefore, we conservatively argue that the 
/bnc fractions obtained in this paper have to be regarded 
as upper limits. 



6 DISCUSSION 

The results presented in the previous section suggest that 
in as much as 25 to 40 percent of all haloes, the bright- 
est galaxy is a satellite rather than a central galaxy. In or- 
der to put these numbers in perspective, we compare them 
to predictions from halo occupation statistics (Section 16. ip 
and from two semi-analytical models of galaxy formation 
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Figure 8. Cumulative distributions of \TZ\ (left panels) and \S\ 
(right panels) of 1000 mock haloes with M = 3 X 10 14 /i _1 M 
populated with galaxies. Black solid curves indicate the distribu- 
tions for haloes with no substructure, while the long-dashed green 
curves, short-dashed blue curves, and dotted red curves indicate 
the distributions for / su b = 1/3, 1/5, and 1/10, respectively. Re- 
sults are shown for the case in which all haloes have substructure 
(upper panels) and 10% of the haloes have substructure (lower 
panels). 



(Section 16. 2|) . We also investigate the impact of non-zero 
/bnc (AO on the halo masses inferred using satellite kine- 
matics (Section 16. 3|) . 



6.1 Comparison with Halo Occupation Statistics 

Using the SDSS r-band luminosity function of Blanton et 
al. (2003), the projected two-point correlation functions of 
Wang et al. (2007), and the galaxy-galaxy lensing data of 
Mandelbaum et al. (2006), C09 constrained the halo occu- 
pation statistics as parameterized by the conditional lumi- 
nosity function (CLF) 



$(L|M) = $con(L|M) + $ sat (L|M) . 



(21) 



which is found to be in good agreement with direct con- 
straints from the Y07 group catalogue (Yang, Mo & van 
den Bosch 2008, 2009a) and with constraints from satellite 
kinematics (More et al. 2009) . We have used this CLF model 
in Section [4.11 to construct our mock group catalogues, ex- 
cept that we imposed that the central galaxy is always the 
BHG; if a satellite luminosity was drawn from $ sat (I/|Af) 
that was brighter than the central luminosity, drawn from 
$cen(i|M), it was rejected. Later, we then switched the lu- 
minosities of the central and that of its brightest satellite in 
a fraction /bnc of all groups. We now use the CLF to 'pre- 
dict' the fraction /bnc as function of halo mass, by making 
the assumption that the luminosities of satellite galaxies are 
independent of the luminosity of the central galaxy in the 
same halo. In that case, the probability that a random satel- 
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lite galaxy in a halo of mass M has a luminosity smaller than 
that of its central is simply given by 



P(L 



sat ^ L c , 

si 



f 

Jo 



|M) = 

;"$ sat (L|M)dL 



(22) 



$ aat (L|M)dL 



$cen(icen|M)di c 



where the integral accounts for the scatter in the relationship 
between central galaxy luminosity and halo mass, and L m i n 
is the minimum luminosity considered. In a halo with N aa ± 
satellites, the probability that the central galaxy is also the 
brightest galaxy is simply [P(L sat < L ccn | M)] Bat . Hence, 
we obtain that 

oo 

/bnc(M) = 1 - P(jV sat |M) [P(L sat < L ccn \M)] N ^ 

JV sat =l 

(23) 

where P(JV sat |M) gives the probability that a halo of mass 
M contains N aa .t satellites, which we take to be a Poisson 
distribution whose mean is given by Eqn. ©. The results 
thus obtained from the CLF of C09 are shown as a solid 
line in Fig. [9] For comparison, the solid triangles with er- 
ror bars are the constraints on /bnc (A/) obtained from the 
SDSS group catalogue as described in the previous section. 
Although the CLF predicts that /bnc increases with increas- 
ing halo mass, in qualitative agreement with the data, the 
values of /bnc (M) inferred from the data are much larger 
than those 'predicted' by the CLF. Note that the CLF pre- 
diction is for all haloes, rather than for haloes (or groups) 
with iV sa t > 3, which is the richness threshold used in our 
analysis of SDSS groups. However, this has little to no im- 
pact. We have verified that our constraints on /bnc do not 
change significantly if we use a different richness threshold. 
Furthermore, from the CLF, the probability of a halo host- 
ing fewer than three galaxies is low for most of the mass 
range considered here: P(N sa .t < 3) > 0.05 only for haloes 
with log M < 12.75. 

There are a number of possible explanations for why 
the CLF prediction is not consistent with the data. First 
of all, we emphasize that the CLF is not designed to 'pre- 
dict' fBNc(M). Most of the statistics used to constrain the 
CLF depend only very weakly (clustering and lensing) or 
not at all (luminosity function) on /bnc(M). Secondly, it 
is still possible that the CLF is correct, but that the addi- 
tional assumption that the satellite luminosity is indepen- 
dent of the luminosity of the central is not correct, that is, 
$ aa t(P S at|M, Lcen) / $ sa t (Lsat | M) . This could come about, 
for example, because of galactic 'cannibalism': those haloes 
in which the central has recently cannibalised a bright satel- 
lite will have an excessively bright central, and are less likely 
to have a bright satellite. It remains to be seen whether 
a model for $ sa t(L S at|M, L ccn ) can be found that yields a 
/bnc (AO in better agreement with the data, and simulta- 
neously obeys the CLF constraint that 

$ sa t(£ 8 at|M) = J $ sa t(isat|M, icon) $cen(£cen | M) di C on . 

(24) 

Of course, it is also possible that the discrepancy re- 
flects an actual failure of the CLF. One possible modifica- 
tion, which will have little impact on the luminosity func- 
tion, clustering, galaxy-galaxy lensing, and mock group cat- 
alogues, is a modification in the shape of $ sa t(L|M) (eqn. [9]) 




Figure 9. Probability that the most luminous satellite galaxy is 
brighter than the central galaxy in a halo, as a function of halo 
mass (Eqn. 1231 1. Result is shown for three slopes in the satel- 
lite CLF (Eqn. EJ: s = 1 (dotted curve), the fiducial s = 2 
(solid curve), and s = 3 (dashed curve). For comparison, the 
/BNc(A^halo) result from Figure [7] (solid triangle points) is also 
shown. The predictions from the MORGANA semi-analytic model 
(open triangles) and Croton et al. (2006) semi-analytic model 
(open squares) are also shown, with Poisson errors. 



at the bright end. As mentioned in Section R.ll C09 adopted 
a functional form for which & ss ,t(L\M) oc exp [— (L/L^ t ) s ] 
at the bright end, with s = 2. The exact value of s, though, is 
poorly constrained by the data, but has a significant impact 
on fBNc(M). This is illustrated by the dotted and dashed 
curves in Fig. [9l which correspond to s = 1 and s = 3, 
respectively. Clearly, decreasing s increases the expectation 
value for the luminosity of the brightest satellite, and hence 
the fraction of haloes for which the central galaxy is not 
the BHG. For s = 1 the CLF 'predicts' a /bnc(M) in good 
agreement with the data, but only for M > 5 x 10 13 /i _1 M©. 
For less massive haloes, the value of /bnc inferred from the 
SDSS group catalogue is still too high compared to the CLF 
prediction, but only by about 2a. 

Another parameter that has a significant impact on 
/bnc(-W) is the ratio Q = L* at /L cen . Motivated by the 
CLF inferred from the SDSS group catalogue by Yang et 
al. (2008), C09 adopted Q = 0.562, with no dependence 
on halo mass. However, Hansen et al. (2009), using the 
maxBCG group catalogue of Koester et al. (2007), found 
that Q depends on halo mass and may be as small as ~ 0.15 
for M ~ 10 14 /i _1 M Q . In order for the CLF to yield /bnc(M) 
in rough agreement with our findings, though, we need a 
value of Q that is larger than 0.562, (i.e., for s = 2 we 
require Q ~ 0.7). These results warrant a more thorough 
investigation into the exact values of s and Q as function of 
halo mass. 

Finally, we have verified that the CLF predictions for 
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fBNc(M) do not depend significantly on our choice for the 
lower luminosity cut-off, I/ m i n - 

6.2 Comparison with Semi- Analytical Models 

We now compare our results to the predictions of two semi- 
analytical models (SAMs) of galaxy formation; the MOR- 
GANA model of Monaco, Fontanot & Taffoni (2007), as up- 
dated by Lo Faro et al. (2009), and the SAM of Croton et 
al. (2006). 

Both SAMs adopt flat ACDM cosmologies, albeit with 
slightly different values for the cosmological parameter^]. 
Although both SAMs include treatments of cooling, star for- 
mation, feedback from supernovae and active galactic nuclei, 
mergers, starbursts and disk instabilities, the actual imple- 
mentations of these physical processes are substantially dif- 
ferent (see the original papers for details). Yet, as shown in 
Fig. [9] (open symbols), they predict fairly similar values for 
/bnc(A^)- In qualitative agreement with the data, both 
SAMs predict that /bnc increases with halo mass. In addi- 
tion, from studying the satellite BHGs in the MORGANA 
model, we find that the majority of the more massive satel- 
lites represent a recently (since z ~ 0.2) accreted population, 
in most cases linked to the last major merger experienced 
by the host halo. In other words, infalling massive satellites 
could be a contributing cause of a nonzero /bnc. 

As with the halo occupation statistics, the predictions of 
both models are significantly lower than the data. In both 
SAMs, satellite galaxies are 'strangulated' after being ac- 
creted by their host galaxies, so that they do not partici- 
pate in the cooling flow of their parent halo. However, it 
has been pointed out that the standard, instantaneous im- 
plementation of this strangulation causes an over-quenching 
of satellite galaxies (e.g., Baldry et al. 2006; Weinmann et 
al. 2006; Kimm et al. 2009). It has been suggested that this 
over-quenching problem can be avoided by adopting a longer 
time-scale for strangulation (e.g., Kang & van den Bosch 
2008; Font et al. 2008; Weinmann et al. 2009). This is likely 
to allow satellite galaxies to continue forming stars for some 
period, which will increase their stellar mass, and thus also 
/bnc - The effect, though, is likely to be small. Another effect 
that may cause the SAMs to underpredict /bnc is the fact 
that they often adopt dynamical friction time-scales that 
are too short (Boylan-Kolchin et al. 2008; Wetzel & White 
2010). This results in (massive) satellites being accreted too 
rapidly, and hence in an underestimate of /bnc- Finally, 
the models don't take proper account of stellar mass strip- 
ping due to tidal forces, which will make satellites less mas- 
sive. As emphasized in a number of recent papers, this tidal 
stripping of satellite galaxies is an important ingredient of 
galaxy formation (Monaco et al. 2006; White et al. 2007; 
Conroy, Ho & White 2007; Conroy, Wechsler & Kravtsov 
2007; Kang & van den Bosch 2008; Yang, Mo & van den 
Bosch 2009b; Pasquali et al. 2010). It remains to be seen 



4 We do not believe that these small differences will have a sig- 
nificant impact on the predictions of /bnc(A^)- 

5 As in our analysis of SDSS galaxy groups, the predictions of 
/bnc(A^) from the SAMs are the fractions of haloes in which a 
satellite galaxy is the most luminous galaxy. The SAMs' predic- 
tions are almost exactly same for the fractions of haloes in which 
a satellite is the most massive galaxy (in terms of stellar mass) . 



to what extent these additions and/or modifications of the 
semi-analytical models impact /bnc(A^). For the moment, 
we conclude that the /bnc(M) inferred from our analysis 
of SDSS galaxy groups is uncomfortably high compared to 
predictions from both halo occupation statistics and from 
semi-analytical models of galaxy formation. 

6.3 Implication for Satellite Kinematics 

Studies that attempt to infer the masses of dark matter 
haloes using the kinematics of satellite galaxies always as- 
sume that the CGP is valid (e.g., Zaritsky et al. 1993; McKay 
et al. 2002; van den Bosch et al. 2004; More et al. 2009). 
Since a typical central only has a few satellites, one nor- 
mally stacks many centrals together in order to obtain suf- 
ficient signal-to-noise to measure a reliable satellite velocity 
dispersion, a sat . With sufficiently large galaxy redshift sur- 
veys, one can measure <r sat as a function of the luminosity (or 
stellar mass) of the central galaxies (i.e., one stacks centrals 
in narrow bins in luminosity or stellar mass). 

Consider a halo with A sa t satellites with velocities u sa t,i 
with respect to halo centre, and let the central be stationary 
at the centre of the halo (i.e., v ccn — 0). The satellite velocity 
dispersion measured with respect to the central is: 



2 

^truc 



1 

A sat 



2 

sat.i ■ 



(25) 



Without loss of generality, assume that satellite number 1 is 
misidentified to be the central. The measured velocity dis- 
persion in this case will be: 



E^ s; 



f sa t,l) +{v - ^sat.l)' 



(26) 



The first term sums over the remaining A sat — 1 true satel- 
lites, while the last term is the contribution from the true 
central. Using that (u sa t,i) = 0, and that (v ccn ) = (v^ cn ) — 0, 
it is easy to show that the velocity dispersion measured from 
a stack of many systems, each with A aa t satellites, is equal 
to 



1 + /bnc (1 - iV s - 



(27) 



where /bnc is the fraction of systems in which the central 
is misidentified. Note that when the number of satellites per 
central is one, the measured velocity dispersion is identical 
to the true one, independent of /bnc, that is, the fact that 
some satellites are misidentified as centrals has no impact 
on the inferred satellite kinematics(f] However, in the limit 
A sat — > co, one has that <7 meas = vT+ /BNCCtruc- Since the 
inferred halo mass M oc erf at , this implies that the mass of 
the halo will be overestimated by a factor (1+/bnc) 3//2 - Note 
that Eqn. (|27|l is strictly valid only for stacks of systems that 
all have the same number of satellites. In reality, however, 
there will be scatter in the number of satellites per host in 
the stack. In what follows we assume that we do not make 
a large error if we adopt Eqn. (127|l but with A sat replaced 
by the average number of satellites per central, (A sa t). 



6 This statement ignores the fact that a non-zero /bnc a l so i m_ 
pacts the interpretation of the stacking parameter (i.e., the lumi- 
nosity of the alleged 'centrals') and may result in more interlopers. 
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Figure 10. Fraction /bnc °f haloes in which the brightest galaxy 
is not the central one, as a function of the luminosity of the bright- 
est halo galaxy. The solid line shows /bnc(^bhg) f° r the 'ran- 
dom' mocks, and the dotted line shows the fraction for the 'condi- 
tional' mocks, in which whether a halo has a satellite brighter than 
the central is conditional on the luminosity of the central galaxy. 
The luminosities correspond to r-band absolute magnitudes, such 
that logL = 9.5 corresponds to M r as —19, and logL = 10.3 to 
M r as —21, etc. 



For relatively faint centrals, {N sa t) — 1, and ir mcas ~ 
Ctrue independent of the value of /bnc- However, at the 
bright end (iV sat ) becomes substantially larger than unity, 
and a non-zero /bnc may cause a significant overestimate 
of the true cr sa t. For example, in their analysis of satellite 
kinematics, More et al. (2009) have (N sa , t ) ~ 10 for their 
brightest host bins. If /bnc for these bins is comparable to 
fBNc(M) at the massive end, (i.e., /bnc ~ 0.4), their in- 
ferred halo masses around bright host galaxies will be over- 
estimated by a factor ~ 1.6. 

However, since the halo mass luminosity relation for 
central galaxies is not one-to-one, it is not trivial to in- 
fer /bnc for a certain bin in host galaxy luminosity, Lhost, 
from /bnc (M)- In fact, one can construct mocks in which 
/BNc(ihost) is very small at the bright end, even when 
fBNc(M) is large at the massive end, by making the criterion 
for a halo having a satellite brighter than the central condi- 
tional on the luminosity of the central galaxy. We illustrate a 
particular case of this in Fig. [10] where we show /bnc(^bhg) 
for two types of mocks that have identical fBNc(M) = 0.4. 
These are constructed as follows. We start by populating 
the dark matter haloes in our lOO/i" 1 Mpc simulation box 
described in Section 0] with galaxies of different luminosities 
using the CLF of C09 (see Section |4T1 for details). If the 
luminosity of a satellite is brighter than that of its central, 
a new luminosity is drawn until it is fainter than that of 
the central, so that /bnc = 0. For the first set of mocks, we 
switch the luminosities of the central and that of its bright- 
est satellite in a random fraction of 40 percent of all haloes, 
and we measure the resulting /bnc(£bhg)- We refer to these 



as the 'random' mocks. The mean and scatter obtained from 
ten realizations are shown as the solid line in Fig. 1101 As ex- 
pected, /bnc(£bhg) — 0.4, independent of Lbhg- For the 
second set of mocks, we again start from the mocks with 
/bnc = 0, but we now switch the luminosities of brightest 
satellite and central in those haloes that meet the criterion 



$ cen (L|M)dL < 0.4, 



(28) 



that is, in those haloes in which the luminosity of the central 
falls in the lower 40 percentile of its distribution, $ cen (£|M). 
We refer to these as the 'conditional' mocks, and similar to 
the 'random' mocks they thus have /BNc(Af) = 0.4. Yet, 
their /bnc(Jjbhg), indicated by the dotted curve in Fig. 1101 
are very different: they decrease from ~ 0.4 at the faint 
end to almost zero at the bright end. Hence, if the proba- 
bility that a central is not the BHG is conditional on the 
luminosity of the central, which seems to be a reasonable 
assumption, then even a relatively large /bnc (M) may have 
negligible impact on the halo masses inferred from satellite 
kinematics. 

We thus conclude that a proper assessment of the im- 
pact of a non-zero /bnc (Af ) on the halo masses inferred from 
satellite kinematics requires an independent assessment of 
the fraction /bnc as function of the luminosities of the host 
galaxies. 



7 CONCLUSIONS 

It is generally assumed that the central galaxy in a dark 
matter halo, that is, the galaxy with the lowest specific po- 
tential energy, is also the brightest halo galaxy (BHG) and 
that it resides at rest at the centre of the dark matter poten- 
tial well. This central galaxy paradigm (CGP) is an essential 
assumption made in various fields of astronomical research 
(e.g., satellite kinematics, gravitational lensing, both weak 
and strong, halo occupation modelling). 

In this paper, we have used a large galaxy group cata- 
logue, constructed from the SDSS DR4 by Yang et al. (2007) , 
in order to test the validity of the CGP. For each group we 
compute two statistics, 1Z and 5, which quantify the off- 
sets of the line-of-sight velocities and projected positions of 
brightest group galaxies relative to the other group mem- 
bers. By comparing the cumulative distributions of \1Z\ and 
|<S| to those obtained from detailed mock group catalogues, 
we have tested the null- hypothesis, Ho, that the CGP is cor- 
rect; hypothesis "Hi, according to which central galaxies are 
BHGs but have a non-zero velocity with respect to the halo 
centre, parameterized by a non-zero velocity bias, fe ve i; and 
hypothesis H2, according to which central galaxies reside at 
rest at the centre of the halo's potential well, but are not 
the BHGs in a fraction /bnc of all haloes. 

In agreement with vdB05, who only used the 1Z statis- 
tic, we show that the null-hypothesis is strongly ruled out. 
However, contrary to vdB05, who argued that the data are 
consistent with a non-zero 6 vc i, we show that Hi fails to si- 
multaneously match the \TZ\ and |<S| distribution^. Rather, 



7 Since vdB05 did not consider the spatial offsets (i.e., the S 
statistic), they were unable to notice this problem for the Hi 
hypothesis. 
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we have shown that the data is consistent with hypothe- 
sis H2, indicating that BHGs are not central galaxies in 
a non-negligible fraction of all haloes; in particular, we 
find that /bnc increases from ~ 0.25 in low mass haloes 
(10 12 /i _1 M < M < 2 x 10 13 /i _1 M s ) to ~ 0.4 in massive 
haloes (M > 5 x 10 13 /i _1 M ). 

Considering combined models, in which BHGs are not 
central galaxies in a non-zero fraction /bnc of all haloes, and 
central galaxies have a non-zero value for their velocity bias 
6 vc i, we find that the data can at most accommodate a small 
amount of velocity bias (6 ve l ~ 0.2). We emphasize, though, 
that a non-zero velocity bias is not required by the data (see 
vdB05 for a discussion of possible physical explanations for 
a non-zero 6 V ei). 

Our main result is that the fraction /bnc is significant 
and increases with halo mass. Since some authors have as- 
sumed that the most luminous (or most massive) galaxy in 
a system is the central galaxy, our result that /bnc is some- 
what large in galaxy groups and larger (43^%) in more 
massive clusters may seem surprising. Nevertheless, our re- 
sults are consistent with other studies. For example, von 
der Linden et al. (2007) find that in 343 of their 625 clus- 
ters (« 55%), the identified BCG is not the 'mean' galaxy 
(which lies at the centre of the cluster's density field). Co- 
ziol et al. (2009) find that in about half of their 452 clusters, 
the BCGs have a median peculiar velocity greater than one 
third of their clusters' velocity dispersion. Finally, in the Lo- 
cal Group, although the Milky Way and Andromeda both 
have their own systems of satellites, they could be considered 
part of a single group, in which Andromeda would be iden- 
tified as the BHG (van den Bergh 1999). The two galaxies 
are expected to eventually merge, although neither is clearly 
the 'central' galaxy. 

In order to put our constraints on /bnc (AO in per- 
spective, we have compared them to predictions from 
semi-analytical models (SAMs) for galaxy formation and 
from halo occupation statistics. Both the SAM of Cro- 
ton et al. (2006) and that of Lo Faro et al. (2009) 
predict 0.1 < /bnc ^ 0.2 in the halo mass range of 
10 13 /i _1 M < M < 10 15 /i _1 M , significantly lower than 
the /bnc (M) inferred here from our SDSS galaxy group cat- 
alogue. 

We can also use the CLF model of C09 to predict 
fBNc(M), if we assume that the luminosity of a satellite 
galaxy is independent of the luminosity of its correspond- 
ing central. Although the C09 halo occupation model ac- 
curately matches the SDSS r-band luminosity function of 
Blanton et al. (2003), the projected two-point correlation 
functions of Wang et al. (2007), and the galaxy-galaxy lens- 
ing data of Mandelbaum et al. (2006), this model also pre- 
dicts a /bnc (AO that is far too low. We have shown that 
one can increase the predicted /bnc(AO using some small 
modifications of the CLF, but it remains to be seen how 
these modifications impact the clustering and lensing pre- 
dictions. For completeness, we emphasize that halo occu- 
pation models based on the abundance-matching method 
(e.g., Vale & Ostriker 2004; Conroy, Wechsler & Kravtsov 
2006, 2007; Shankar et al. 2006; Guo et al. 2010; Moster 
et al. 2010) 'predict' that /bnc = 0, by construction, un- 
less they assume a non-zero scatter in the relation between 
central galaxy luminosity and halo mass. Typically a larger 
amount of scatter will imply a larger /bnc- It remains to 



be seen whether the amount of scatter required to match 
the /bnc (A/) inferred here is consistent with independent 
constraints, such as those obtained from satellite kinematics 
(More et al. 2009). All in all, we conclude that the con- 
straints on /bnc (AO obtained in this paper are uncomfort- 
ably high compared to predictions from galaxy formation 
models and halo occupation statistics. One possible expla- 
nation may be that dark matter haloes have substructure, 
something that we have ignored in our analysis. Although 
simple tests suggest that the correlations in the phase-space 
parameters of satellite galaxies due to substructure are not 
strong enough to significantly impact our conclusions, we 
caution that more detailed studies are required to confirm 
this. 

Our results have important implications for various ar- 
eas in astrophysics. In particular, we have shown that a non- 
zero /bnc may cause an overestimate of halo masses inferred 
from satellite kinematics. Although the effect is expected to 
be negligible for faint host galaxies, because they typically 
only have of the order of one satellite per host, the overes- 
timate can be significant at the bright end. The exact im- 
pact, though, depends on the fraction /bnc as function of 
the host luminosity, which may be very different from that 
as a function of halo mass. For example, if BHGs are not 
centrals in the fraction /bnc (AO of haloes of mass M that 
host the faintest centrals, then /bnc(L) drops towards zero 
at the bright end, and the impact on satellite kinematics is 
negligible. 

Additional methods and analyses that may be affected 
by a non-zero /bnc are the following: 

• The inference of halo masses from weak gravitational 
lensing. Similar to satellite kinematics, weak lensing studies 
often rely on stacking the lensing signals of many clusters 
and groups binned by mass-correlated observables such as 
richness and luminosity. When interpreting such data, it is 
generally assumed that BHGs coincide with the centres of 
the dark matter haloes (e.g., Mandelbaum et al. 2006; John- 
ston et al. 2007; Sheldon et al. 2009a,b; Corless & King 2009; 
C09). Since satellite galaxies yield a different lensing signal 
than centrals, at least on small to intermediate scales (see 
e.g., Yang et al. 2006), a non-zero /bnc may have a signif- 
icant impact on the lensing signal (Johnston et al. 2007), 
by resulting in underestimates of the richnesses and lensing 
profiles, for example. C09 used the CLF described in Sec- 
tion 14.11 to model the galaxy-galaxy lensing signal of Man- 
delbaum et al. (2006). As shown in Section 16.11 this CLF 
predicts /bnc (AO significantly lower than inferred here. It 
remains to be seen to what extent this may effect their in- 
terpretation of the lensing data. 

• Analyses of the power spectrum of luminous red galax- 
ies (Tegmark et al. 2006; Reid et al. 2010), which can be 
used to constrain cosmological parameters. In tests with 
mock galaxy catalogues, Reid et al. (2010) find that a 
fraction /bnc of 0.2-0.4 results in an angle- averaged LRG 
power spectrum that is damped by 2 — 4% on scales of 
k ~ 0.1/iMpc _1 , and the effect is larger at larger k. However, 
when analyzing the observed power spectrum with these 
mock catalogues, the effect of /bnc on the recovered cos- 
mological parameters is relatively small. 

• Measurements of the radial number density distribution 
of satellite galaxies, n sa x(r\M), in haloes of mass M. Sev- 
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eral studies that have measured n aat (r|M) using groups and 
clusters have assumed that the halo centre coincides with 
the location of the BHG (e.g., Carlberg, Yee & Ellingson 
1997; Collister & Lahav 2005; Hansen et al. 2005; Yang et 
al. 2005b) . If, instead, the BHG is a satellite galaxy, this will 
result in an underestimate of the concentration of n sat (r-). 
Hence, a non-zero /bnc(A^) may cast doubt on the claim, 
made by several of these studies, that satellite galaxies are 
less centrally concentrated than the dark matter. On the 
other hand, Lin, Mohr & Stanford (2004) used the centre 
of the X-ray emission as the centre of the dark matter halo, 
rather than the location of the BHG, and came to a similar 
conclusion. 

• Comparisons of the properties of central and satellite 
galaxies. A number of recent studies have used galaxy group 
catalogues to split the galaxy population into centrals and 
satellites, and to compare their properties (e.g., Weinmann 
et al. 2006, 2009; Skibba et al. 2007; van den Bosch et 
al. 2008; Pasquali et al. 2009, 2010; Hansen et al. 2009; 
Kimm et al. 2009; Skibba 2009; Guo et al. 2009). Since all 
of these studies have assumed the brightest group galaxy to 
be the central galaxy, they are likely to have underestimated 
the true differences (i.e., a non-zero /bnc blurs the actual 
differences between centrals and satellites). 
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